Anderson disorder in graphene nanoribbons: A local distribution approach 
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Disorder effects strongly influence the transport properties of graphene based nanodevices even 
to the point of Anderson localization. Focusing on the local density of states and its distribution 
function, we analyze the localization properties of actual size graphene nanoribbons. In particular 
we determine the time evolution and localization length of the single particle wave function in 
dependence on the ribbon extension and edge geometry, as well as on the disorder type and strength. 
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I. INTRODUCTION 

Disorder effects in graphene are of particular impor- 
tance on the account of its two-dimensional (2D) lattice 
structure. The single parameter scaling theory predicts 
that in 2D systems arbitrary weak disorder leads to An- 
derson localization of the single particle wave function 1 . 
For graphene it has been argued that due to the lin- 
ear dispersion in the vicinity of the band center the one 
parameter scaling theory does not hold. The problem 
of Anderson localization in graphene is therefore heavily 
debatedlM. 

Accessing Anderson localization both theoretically and 
experimentally, the local density of states (LDOS) is a 
central quantity. By means of the local distribution ap- 
proach, the distribution of the LDOS may be used to 
distinguish localized from extended states^i Nowa- 
days, the LDOS can be directly measured by scanning 
tunneling spectroscopy experiments ^ 1Q i 11 . 

An ordered, infinite graphene sheet is a zero-gap 
semiconductor with a linear density of states near the 
charge neutrality point 12 . Cutting a graphene nanorib- 
bon (GNR) of finite width out of such a sheet, additional 
aspects have to be considered. First, the finite number of 
transverse atoms causes quantum confinement, where the 
presence of the edges leads to a symmetry breaking. Sec- 
ond, lattice defects or targeted implementations of e.g. 
boron (B7) clusters 13 result in random samples. Thereby 
the range of the disorder is of great importance 14 : For 
long-range disorder, as caused by ripples in the graphene 
sheet, the two independent corners of the Brillouin zone 
are untangled and long-wavelength excitations can be 
modeled by an effective Dirac equation. If the scattering 
potential is short ranged, however, inter-valley scatter- 
ing between the two inequivalent Dirac cones becomes 
possible. Third, the finite extension (aspect ratio) of the 
GNR introduces a new length scale being absent in infi- 
nite graphene sheets. Actually we expect metallic behav- 
ior of disordered quasi- ID GNRs if the localization length 
becomes comparable or even larger than the longitudinal 
ribbon siz o 15 i 16 i 17 . 

To address these questions, in this work we investigate 
the electronic structure and the localization properties of 



disordered GNRs by means of unbiased numerical tech- 
niques. Thereby we focus on the interplay of disorder, 
boundaries effects and GNR geometry. Particular as- 
pects of various kinds of disorder in GNRs have been in- 
vestigated previously in the literature 15 1 16 1 17 1 18 1 19 1 20 i 21 1 22 . 



II. MODEL AND METHODS 



To this end we consider the tight-binding Hamiltonian 
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on a honeycomb lattice with N sites, including hopping 
between nearest neighbors (ij) only. Drawing the on-site 
potentials e$ from the box distribution 



pU = -0(7/2-1* 
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we introduce (short ranged) Anderson disorder 23 . We 
distinguish between bulk (75) and edge (j e ) disorder, 
when all on-site potentials are subjected to p[ci] or only 
those at the edge sites. We consider quasi- ID GNRs of 
finite widths with open (periodic) boundary conditions 
in transverse (longitudinal) direction. Depending on the 
orientation of the GNRs with respect to the honeycomb 
lattice, zigzag or armchair geometries will be realized (see 
panels on top of Fig. [1]). 

The local properties of site i of a sample with broken 
translational invariance are reflected in the LDOS, 
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Recording the probability density function f[pi] for many 
different sites {i} of a given sample and different sample 
realizations {e^} restores translational invariance on the 
level of distributions: The shape of f[pi] is determined 
by p[ei] but independent of {i} and {q} 7 . For extended 
states, f[pi) is strongly peaked around the mean DOS, 



= {pi) 



(4) 



2 



independent of the system size, whereas for localized 
states f[pi] exhibits a log- normal distribution that be- 
comes singular for increasing system sizes 24 . Normalizing 
the LDOS to p me allows for a detection of the localiza- 
tion properties by performing a finite size scaling for the 
LDOS distribution. More conveniently, the typical DOS 



Pty 



(5) 



monitors the changes in the LDOS distribution. While 
for TV — > oo an extended state is characterized by finite 
values of p me and p ty , for localized states p me is finite but 

Alternatively, the recurrence probability Pr(£ — > oo) 
also reveals the localization properties of the system 25 . 
While in the thermodynamic limit Pr ~ \/N — > for 
extended states, localized states are characterized by a 
finite value of Pr. Starting from a localized wave packet, 
we are able to calculate the time dependent local particle 
density, 

I N 2 

m(t) = m^t)\ 2 = I e-^imim)^™) , (6) 

rn=l 

by expanding the time evolution operator into a finite 
series of Chebyshev polynomials 26,27 . The above local 
distribution approach also applies to rti(t). But since an 
initial state in general contains contributions of the whole 
spectrum, examining rii(t) does not allow for an energy 
resolved investigation of localization as by the LDOS. 
Instead it provides a tool for a global examination of 
the spectrum with relevance for possible measurements. 
Note that a finite overlap of just one extended state with 
the initial state leads to a complete spreading of this state 
after some time. 



III. NUMERICAL RESULTS 

A. Local density of states 

Compared to the band structure of an infinite 2D 
graphene sheet, the DOS of finite GNRs is character- 
ized by a multitude of van Hove singularities (see top 
panels of Fig. [1]). For zigzag GNRs, the strong signa- 
ture at E = indicates the high degeneracy of the edge 
states 28 . In contrast, armchair GNRs with N a — 3n or 
N a = 3n + 1 are gapped around E = 0. This finite size 
gap tends to zero as N a — > oo. The resulting metallicity 
for N a = 3n + 2 is an artifact of the NN tight binding 
approximation, however, and vanishes if next- and third 
NN are taken into account^. For other values of N a a 
longer ranged hopping slightly modifies the gap size but 
does not change the fundamental behavior. Note that 
even for vanishing Anderson disorder the LDOS varies 
for different bulk sites according to their relative posi- 
tion to the ribbon edges. Symmetry considerations show 
that there are N z (N a /2) inequivalent lattice sites in or- 
dered zigzag and armchair GNRs. Therefore mean and 




FIG. 1: (Color online) Mean (solid red) and typical (dashed 
blue) DOS for zigzag (left column, N z = 6) and armchair 
(right column, N a = 10) GNRs of width W = 1.1 nm. Top 
panels: Ordered case. Lower panels: In each (2 x 2) block, 
we compare for a fixed value of disorder the influence of bulk 
disorder (75, upper rows) to edge disorder (7^ bottom rows). 
To illustrate the localization properties, in each panel pt y is 
given for L = 213 (1064) nm by lightblue dot-dashed (dark- 
blue dashed) lines. These system sizes correspond to 10000 
(50000) lattice sites for the armchair and 10392 (51960) for 
the zigzag case. Disorder averaging was performed over 10 5 
realizations. In the longitudinal direction periodic boundary 
conditions (pbc) are applied. 



typical DOS do not coincide even for 75 = (see, e.g., 
the band center of the zigzag GNR). 

If disorder comes in, localized states emerge in the 
band gap of the armchair GNRs, and above a critical 
disorder strength the gap is filled completely. The local- 
ization properties of the states can readily be seen from 
the system size dependence of p ty . A tendency towards 
reduced values of p ty for increasing system sizes indicates 
localization for both GNR geometries and all energies. 
While this localization effect arises for bulk disorder al- 
ready at 75 jt = 2, an edge disorder strength of j e /i = 2 
is still too weak to localize the wave function on GNRs of 
L = 213 nm size as indicated by the approximate equality 
of pty and p me . A substantial reduction of p ty is only ob- 
served for larger systems (L = 1064 nm) which indicates 
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FIG. 2: (Color online) Gap size A a for armchair GNRs as a 
function of ribbon width N a = 3n + 1 for edge disorder (main 
panel) and bulk disorder (inset). Results are based on the 
averaged DOS for GNRs with a length of 1000 atoms using 
4096 realizations of disorder. 



localization on a larger length scale. Obviously, zigzag 
GNRs are less sensible to edge disorder than armchair 
GNRs since this geometry has only half the number of 
(disordered) edge sites. The different edge geometries are 
only of importance if the disorder is weak. For strong 
disorder, jb/i = 4, the results for armchair and zigzag 
GNRs coincide almost exactly. 

As stressed above, there are three branches of gap sizes 
depending on mod(7V a ,3). In Fig. [2] we focus on N a = 
3n + 1 and examine the influence of both bulk and edge 
disorder on the gap size A a in dependence on the ribbon 
width. For our finite system we calculate A a as 



A a /2 



Pme (E)dE=±. 
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A finite-size analysis shows that upon increasing the rib- 
bon width the gap narrows for any bulk disorder 75. In 
contrast, for edge disorder we observe a non-monotonic 
behavior that can be explained by the competition of 
two effects: Increasing the width of the GNR on the one 
hand weakens the influence of the disorder as the ratio 
of edge to bulk sites decreases. An increasing number of 
lattice sites, on the other hand, reduces the finite size ef- 
fects and closes the gap. Thus, for j e > 2f, the gap first 
broadens when the GNR width is increased, and then 
converges to the gap size of the ordered system, which 
finally vanishes in the limit N a — > 00. Similar studies 
for a different type of edge disorder, in which sites are 
randomly removed from the ribbon edges, can be found 
in the literaturei^ 2 ^. 

To get further insight into the nature of the eigenstates 
of GNRs and substantiate our conclusions about their lo- 
calization properties, we show the LDOS in the band cen- 
ter in Fig. [3j The magnifying inset for the ordered case 
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FIG. 3: (Color online) Normalized LDOS at the band cen- 
ter (pi/ p m e)\E=o for particular zigzag GNRs. In addition 
to contrasting bulk and edge disorder, we compare in the 
left (right) column the influence of the aspect ratios L x W = 
31.4x13.5 (31.4x3.3) nm 2 , corresponding to 256x64 (256x16) 
sites. Results obtained by exact diagonalization. 



shows the alternating structure of the edge states which 
are distinctive for the band center of zigzag GNRs 28 . 
In the presence of weak edge disorder, the checkerboard 
structure of the amplitudes persists in the bulk, while 
near the edges regions with significantly enhanced ampli- 
tudes emerge. The A-B sublattice structure is no longer 
present for larger j e as can be seen in the lower inset of 
Fig. [3j Here, the sites with vanishing amplitudes form a 
filamentary network in the bulk, caused by the influence 
of the disordered edges. For bulk disorder, localization 
arises first near the edges of the system in the case of weak 
disorder, while localized states in the bulk of the GNR 
occur only for strong disorder. Varying the aspect ratio 
of the GNRs (right column of Fig. [3]), we may tune the 
relative importance of the edges in the system. Although 
this effect is most pronounced for edge disorder, we ob- 
serve also for bulk disorder such a "renormalization" of 
the disorder strength: A given 75^ causes stronger local- 
ization for narrow GNRs. 



B. Time evolution of the wavefunction 

Figure [H shows the time evolution of an initially lo- 
calized state, as calculated by the Chebyshev method 27 . 
The dynamics of the initial wave packet is characterized 
by a fast spreading process (t < 10 3 to), after which its 
extension does not change anymore, even for very long 
times. Clearly, on individual sites the amplitudes fluc- 
tuate in time, but the overall nature of the state for 
t = IOHq is quasistationary. The localization proper- 
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FIG. 4: (Color online) Time evolution of the normalized par- 
ticle density N\ip(ri)\ 2 on disordered GRNs with zigzag and 
armchair geometries for different values of bulk disorder 75. 
Device dimensions: (1.1 x 213) nm 2 corresponding to 6 x 1732 
atoms (zigzag) and 10 x 1000 atoms (armchair). Times are 
measured in units of the inverse hopping element to = 1/t. 
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is no problem to determine A from an exponential fit 



(r,)| 2 = kHr )| 2 exp 



A 



(8) 



for a given initial state and disorder realization at any 
fixed time. But the such-determined A strongly fluctu- 
ates, both in time and as a function of the chosen ini- 
tial state and disorder realization. The temporal fluctu- 
ations of about 5-10% can be eliminated by time averag- 
ing. Varying the initial state and/or comparing different 
disorder realizations, leads to additional uncertainties of 
about 10-20%. Therefore we show in Fig. [5] sample aver- 
ages over several combinations of initial states and disor- 
der realizations. 

Figure [5] indicates that the influence of the boundary 
(armchair/zigzag) is only of minor importance for the 
localization length. But we observe a pronounced differ- 
ence between bulk and edge disorder, with A > L also 
for large values of j e for most ribbon widths. For any 
fixed disorder strength, a decreasing width of the GNR 
systematically reduces A since the influence of the lateral 
dimension is weakened and the system approaches the 
ID limit. Values of A which are significantly larger than 
half the system size (blue solid line) have to be taken 
with care since a reliable determination of the localiza- 
tion length requires A < L. Clearly, the precise value of A 
in those cases is of minor importance due to the metallic 
behavior of such finite GNRs. A quantitative compar- 
ison of the obtained localisation lengths with estimates 
based on other method o 15 i 16 i 17 suffers from the different 
investigated disorder models. Nevertheless, the orders of 
magnitude match and the general tendencies are repro- 
duced: the impact of disorder increases with decreasing 
ribbon width and the boundary type does not influence 
the localization length significantly for strong disorder. 
The pronounced dependence of the localization length 
on the ribbon type (armchair or zigzag) for the weakly 
disordered case reported in Refs. fl5lfl6l is absent in our 
data. We attribute this to the different disorder models 
used. 



FIG. 5: (Color online) Localization length in dependence on 
bulk (7b) and edge (pf e ) disorder strength for armchair and 
zigzag GNRs. The values are sample averages obtained for 
10 GNRs of L — 213 nm when the state has become quasis- 
tationary. 



ties depend on both disorder strength and edge geome- 
try. Obviously, armchair GNRs are more susceptible to 
the presence of disorder than those of zigzag type. For 
the shown GNRs of moderate length and weak disorder 
(76, e A = 0.5) the localization length is larger than the 
system size and thus the GNR is "metallic" . 

The extraction and quantitative discussion of the lo- 
calization length in narrow GNRs is challenging. There 



IV. SUMMARY 

To conclude, Anderson localization takes place in dis- 
ordered quasi- ID graphene nanoribbons, but taking into 
account the actual device dimensions GNRs can be con- 
ducting at weak disorder strengths. This has been proven 
by calculating the localization length and time evolution 
of single particle states. Within the local distribution ap- 
proach Anderson localization is identified by a log-normal 
distribution of the LDOS that shifts towards zero for in- 
creasing system size. The LDOS is directly measurable 
by scanning tunneling spectroscopy and therefore allows 
for a direct comparison of theory and experiment. 
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